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Some  integrals  are  presented  that  can  be  expressed  in  terms  of  the  Q-function, 
which  is  defined  as 

lo(ax)  , 

and  where  Iq  is  the  modified  Bessel  function  of  order  zero.  Also,  integrals  of  the 
Q-function  are  evaluated.  Some  of  the  integrals  are  generalizations  of  earlier  results, 
but  others  are  new;  all  derivations  are  included.  Extensions  to  related  integrals  are  also 
presented. 
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ABSTRACT 


Some  integrals  are  presented  that  can  be  expressed  in  terms 
of  the  Q-function,  which  is  defined  as 


Q  (a,  b) 


dx  X  exp  l- 


2  2 
X  +  a 


Iq  (ax) 


and  where  Iq  is  the  modified  Bessel  function  oforderzero.  Also, 
integrals  of  the  Q-function  are  evaluated.  Some  of  the  integrals 
are  generalizations  of  earlier  results,  but  others  are  new;  all  de¬ 
rivations  are  included.  Extensions  to  related  integrals  are  als 
presented. 
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SOME  INTEGRALS  INVOLVING  THE  Q-FUNCTION 


INTRODUCTION 

The  performance  analysis  of  phase-incoherent  receivers  in  fading  or 
nonfading  media  requires  evaluating  the  Q  -function.  It  is  defined  as  [Ref.  1, 
Eq.  (16)] 


f  /  2  ,,  2\ 

Q(a.b)  =  j  dx  X  expf-^  ^  jl^(ax) ,  (1) 

where  Iq  is  the  modified  Bessel  function  of  order  zero.  Physically,  the  Q- 
function  gives  the  cumulative  distribution  function  of  the  envelope  of  the  sum  of 
a  sine  wave  and  a  narrowband  Gaussian  process  [Ref.  2]. 

Much  past  work  on  performance  analysis  of  receivers  in  fading  media  hao 
required  evaluating  the  Q-function  or  its  integrals  [Ref.  3-10).  It  is  the  pur  ¬ 
pose  of  this  report  to  compile  past  results  and  generalize  and  augment  them, 
where  possible.  (The  method  for  deriving  every  relation  is  included  (as  an 
appendix)  so  that  the  reader  can  formulate  his  own  generalizations ,  where 
appropriate.)  The  availability  of  closed -form  expressions  for  these  integrals 
(in  terms  of  the  Q-function)  greatly  facilitates  numerical  evaluation,  since 
programs  for  the  Q-function  are  available  [Ref.  11] . 

Some  generalizations  and  extensions  of  the  Q-function  are  also  included. 


GENERAL  RELATIONS  FOR  THE  Q-FUNCTION 

Q(a,  0)  =  1,  Q(a, «»)  =0  (a  <  •«) ,  * 


Q(0,b)  =  exp(  >  /2),  Q(<»,b)  =1  (b  <«*) . 


(2) 


*Derivations  of  these  and  succeeding  relations  are  presented  m  the  ap¬ 
pendix  . 


1 


I 
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IISTEGRALS 

In  the  mtegrals  listed  below,  tjxe  parameters  are  presumed  to  be  real  and 
positive.  However,  the  results  may  be  generalized  to  negative  or  complex 
parameter  values  in  many  cases  by  symmetry  or  analytic  continuation.  For  ex¬ 
ample,  (9)  holds  for  all  complex  a  and  b  and  complex  p,  provided  Re(p^)>0. 
This  follows  because  the  integrand  of  (9^  is  analytic,  aiid  the  integral  is  uni¬ 
formly  convergent  for  Re(p^)  >  0;  see  Ref.  12,  pp.  99-iCO.  The  right-hand 
side  of  (9)  is  also  analytic  for  all  a  and  b  and  for  p  ^  0.  As  an  exao'ple  s-f 
the  use  of  symmetry,  consider  (77),  which  is  derived  for  b  >  0.  Since  is 
an  odd  functiori,  we  can  express  J^(bx)  =  sgn(b)  J]^((b|x)  for  all  real  b.  and 
utilize  the  given  result  by  substituting  |b|  for  b  and  multiplying  by  Kgn(b). 
The  integrals  have  been  checked  numerically. 


BESSEL  FUNCTIONS,  EXPONENTIALS,  AND  POWERS 


J  dx  X  exp(-p^x^/2)  I^lax)  =iexp^^  Q(a/p,bp) 

b  /  2  \ 

f  dx  X  exp(-p^x^/2)  I^(ax)  “■^exp/- J  [1  -  Q(a/p,bp)] 

J  \  n  \2p  ' 

b  /  2\ 

f  2  2  1  /  \ 

1  dx  X  exp(p  X  /2)  I^(ax)  =  — expl - [Q(ia/p,ibp)  -  1] 

n  p  \  2p  / 


n 

^  dx  X 


t/  2  2\  .  /  2 

jxp^'lJL.j  1^(1)  p^)  _  exp|-  ^ 


^  dx  exp(-p^x^/2)  Ij^(ax)  =  [l  -  Q(bp,a/p)] 

b 

(  dx  exp(-p^x^/2)  I  'ax)  =  -■  [exp^^-r:^Q(bD,a/D)  -  ll 

J  ’  H  J 


3 


TR  4297 


b  2 

j  dx  exp{pV/2)  I^(ax)  =  ^  pp^-^jQ(ibp,ia/p)  -  1 

b 

I  dx  exp(pV/2)  I^(bp\)  =-ij[^iexp{-bV/2) 


+  ^exp(bV/2)  iQ(bV)  -  l] 


(X>  r-  O 

C  2  2*^  li  /^\ 

J  dxx  exp(-p  xV2)I^(ax)=~|^aexp^jQ(a/p,bp) 

+  bp^  exp(-b^p^/2)  Ij^(ab) 


j  dx  X  exp(-pV/2)  I^(ax)  =  ~  a  exf/^  jl  -  Q(a/p,bp)| 


2  2  2 
-  bp  exp(-b  p  /2)  I 


^(ab) 


w  ^  2 

j  dx  exp(-pV/2)  Ij,(ax)  =  i  (a^  +  2p^  exp/^i^]  Q(a/p,  bp) 


+-^  exp(-bV/2)  [a  I^(ab)  +  bp^  I^(ab)J  (19) 

b  ^ 

J  dx  X  exp(-p  X  /2)  I^(ax)  =-Y(a^  +  2p^)  exp^\  [l  Q(a/p,bp)l 
0  P  \2p  / 

-  -^exp(-bV/2)  [a  I^(ab)  +  bp^  I^(ab)j  (20) 
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(X) 

j  dx  exp(-pV/2»  jQ(ax)  J^(bx)  =  ^  ll  -  Q(a/p,b/p)l 
0 


<K> 

j  dx  exp{-pV/2)  iQ(ax)  I^{bx)  =  ^  lQ(ia/p,ib/p)  -  1] 
0 


tX) 

j"  di.  exp(-p^x^/2)  I^(ax)  Ij(ax) 
0 


(21) 


(22) 


(23) 


In  (U)  -  (34). 


u=^a(p-s),  and  v  =^a(p+sy. 


(K) 

I  dx  exp(-px)  I^(bx)  =  j  [2Q(u,v)-  exp(-pa)  l^(ab)]  ,  p  >  b 


a 

a 


a 

I  d.\  exp(-HA/  .^(px)  =  a  exp(~pa)  Il^j(pa)  +  Ij(pa)] 


(24) 


J  dx  exp(-px)  iQ(bx)  =  ^  [1  +  exp(-pa)  I^(ab)  -  2Q(u,v)]  ,  p4=  b  (25) 


(2G) 


•  1/ 

I  dx  exp(-px)  =~  [2pQ(u,v)  -  (p+s)  ex.p(-po)  I^(ab)l ,  p  >  b  (27) 


a 

a 


j  dx  exp(-px)  Ij^(bx)  ^'-i—  IP  "  s  +  (p+s)  oxp(-pa)  I^(ab)  -2pQ(u,v)j ,  p  #:b  (28) 


a 

j  dx  cxp(-px)  I^(1J\)  --  ^  [exp(-pa)  |^l■^pa)  Iy(pa)  +  pa  I^(i)a)|  -  1] 


(29) 
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2  2  1 

dx  X  exp(-p  X  /2)  Q{ax,b)  =  |exp|- 


-  exp 


(37) 


1 


2  2  1 
dx  X  exp(p  X  /2)  Q(px,b)  exp 


[1  -  Q(b,pc)] 


(39) 


r 

J 


2  2  1 
dx  X  exp(-p  X  /2)  Q(b,ax)  = 


Q(b,ac) 
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<w  r  2  /  ,  2  2  \ 

f  dx  X  exp(pV/2)  Q(b,ax)  =-^ 

J  P  La  -p  \2(a  -p  )/  ' 

/2  2\  ] 

-  exp(^~)  Q(b,ac)  ,  p<a 


^22’ 


c /a  -p 


r2/22\/ 

f  dx  X  exp(p^x^/2)  Q(b,ax)  =  ^  -j—exp|  2")r~^f" 

5  P  La  -p  \2(a  -p  )/( 

+  exp^~^  Q(b,  ac)  -  1  ,  p  ^ 


1  2’ 


c /a  -p 


C 

f  dxxexp(p^x^/2) Q(b,px)  =  -^Ipc  exp(-b^/2)  Ij^(bpc)  +  b  exp(p^c^/2)  Q(b,pc)  - b] 
J  bo  ... 


Q- FUNCTION,  BESSEL  FUNCTIONS,  EXPONENTIALS,  AND  POWERS 


/  2\  / 

J  dx  X  exp(-pV/2)  I^(GX)  Q(ax,b)  = -Lexpp-^jQ^^— ^ 


(2  2 ’  nrz 

p  +a  Jp  +a 


<V  /  /  y 

(  22y  1  /c\/bp  ac 

I  dx  x  exp(-p  X  /2)  I^(cx)  Q(b,ax)  =—  exrf — ^  Q/  — 7=f— j 

J  p  L  \2p  /  \Jp  +a  p^/p  +a‘ 


2  /  2  ^2  2\  /  , 
a  / c  -b  p  \ ,  /  abc  \ 

2  2  2  2^ )  ^o(^~2) 

p  +a  \2(p  +a  )/  \p  +a/j 


/  2  \ 

j  dx  X  Ij,(ox)  Q(b,ax)  =  -i-  [“‘i-r)  *  ““i©] 

^  a  c  v^a  / 


r  /« 2  \  /  V  -1 

\  dx  exp(-p^x^/2)  I^(bx)  Q(c,ax)  =  ^  expj— ^Jqi  - — 7|=y)"  ^ 

L  \2p  /  \  /p  +a  p  /p  +a  /  _ 
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r  2  2  if  /b^\  /  ab 

I  dx  exp(-p  X  /2)  I  (bx)  Q(ax,c)  =  -  expl  — 

J  L  \2p/  \pjp 


+a  ^p  +a 


/  c  \J  ib  iac 


«>  r  /  2  2\  /  2 

f  2  2,  ^  ,  1  /a  c  \^/  a  c _ pc 

1  dx  exp(-p  X  /2)  I  (acx)  Q(ax,c)  =  ~  expl  — ^IQI  i-g — g-,  ™ ^ 

J  L  \2p  /  \pjp  +a  JP 


/  2  2  2\  /  2  2\  ,  /  ii\| 

\  p  +a  /  vp  +a  /  •* 


w  r  /  ^  \  be  1 

f  dx  Ij^(cx)  Q(b,  ax)  =  ^  ®’‘p('^)  ^o(r)  ■  ^ 

L  \2d.  /  J 

0 

^  r  /  2  \  / 

r  9  2  2.  II  /  C  \^/  &C 

I  dx  X  exp(-p  X  /2)I  (cx)  Q(ax,b)  =—  cexpl — ^IQI  .  ^  ^ 

J  ^  n  \2d  /  \d  /d  +a 


2p  '  'pjp  +a  /p  +a 


,2  /  2  .  2  2x  /  ,.  V 

I  expr'^PU  f  ^  ] 

^  2  2®  PI  2  2  I  ll  2  2 

p  +  a  \2(p  +a  )/  ^  +a  / 


r  /  2  \  , 

f  2  2  2,  1  /c  \^/  bp  _ w 

j  dx  X  exp(-p  X  /2)I  (cx)Q(b,ax)  =“  cexpl  2pl  Fz 

J  n  \2p  /  \  ^p  +a  pjp  +a 


(p  +a  )  \2(p  +a  )'{  \p  +a/  ^  ^^/). 


/  2\ 

(  dx  x^  I  (cx)  Q  (b  ,ax)  -  exp  / Fc  (c^  +a^l/ )  Ijj(y) + 2ab{c^  -a^ )  I  j(^)]  (54) 

i  a  c  \2a  /  ^ 
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Q-function  of  two  linear  arguments,  exponentials,  and  X 

in  (5S,  and  (56,.  a„n  r  = 


m 

I  dx  X  exp(-pV/2)  Q(ax.bx)  :=-L  j^exp^-^^j 


dx  X  exp(-p^x^/ 


p-L  V  2 

%  Q(ax.bx,  .  l)j^,  ,  exp(-i|f),^,abc-",] 

/  2  2a 

-exp(-£^jQ(ac.bc)  -~q(c/^. 


(55) 


(56) 


2  2 
<xQ(ax.bx)=|.Q(ac,bc)+— ^ 

«-» *  2  , 


2(a2-b2,  K(nbc^)+aZj(abo^] 


■  i?-b4v, 

-  2Q(c  min(a,b),  c  max(a,b))j,  a  ^  b 
f  2  r 

J  clx  X  Q(ax,ax)  =  ^  [l  *  eapf-a^o^,  {.^(a^c^  .I^,a^e^}] 


2 

0  ^ 


(57) 


(58) 
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(59) 


(GO) 


(Gl) 


(G2) 


(GG) 
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ERROR  FUNCTION,  BESSEL  FUNCTIONS,  EXPONENTIALS,  AND  x 
In  this  report,  the  error  function  is  defined  as 

X 

4>{x)  =  j  dy(27r)"^/^  exp(-yV2) .  ^ 

w 


w 

J  dx  exp(-p^x^/2)  I^(bx)  <|)(ax)  = 


(67) 
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BESSEL  FUNCTIONS  OF  MORE  COMPLICATED  ARGUMENTS 
EXPONENTIALS,  AND  X 

In  (71)-(80),  s  =  t=  •Jq-c,  and  u= 


13 
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_qb^  \j  /  cb^ 


4(q^-cV  ®\4(q^-c^)/ 


I  dx  exp(-qx^)  jj^(bx)  =  i  j^l ,.  exp/- 

on 

I  dx  exp(-ox^  Jj(bx)  ycx^  =  i  [24^  -  1] 


q  >c 


V 

j  <^?‘=»«=<exp(-rrb~)=^[2uQ(i,^,  ijsjraj) 


b<l 


ir 

j  dx  i  cos  X  exp(-  ~  A  =  2a  [*(- JJ)  -  (2xa)'^/^ 


exp  (-a/2)] 


(V 

j  dx  exp(-qx^)  J^(bx)  I^(cx^)  =  ^  |^exp^- 


4(q-e“) 


c  )/ 


q>  c 


m 


I  dx  oxp(-ox^)  I^,cx^  .  exp(-  g)  -  b  4,(. 
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EXTENSIONS 

We  define  a  generalization  of  the  Q-function  as 

b 

Then  Qj^q  the  standard  Q-function.  Integration  by  parts  of  (86)  two  different 
ways  yields  the  relations  (see  (A-49)  and  (A-50)) 


15 
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APPLICATION 


The  error  probability  fo  transmission  of  M-ary  equicorrelated  signals 
over  a  phase -incoherent  Rayleigh  fading  channel  has  been  considered  by  Shein 
(Ref.  9].  If  we  generalize  to  the  case  where,  in  addition,  a  threshold  must  be 
exceeded  (Ref.  5,  sect.  6],  the  probability  of  correct  detection  is  given  by 

'’m'TT?  j  duuexp(-^^)  jdvvexp(-vV2) 

r/Ji^  0 


Iq(s/X  uv)  (1  -  Q(\/X v,u)] 


M-1 


(96) 


(This  specializes  to  Shein's  result  for  the  threshold  equal  to  zero:;  F  =  0,)  For 
M  =  2,  we  use  (45)  to  evaluate  the  inner  integral  in  (96)  and  obtain 


P 


2 


1  -X 
1  +/J 


du  u  exp 


r/v/Tx 


Ql 


Xu  u  \ 


(97) 


This  integral  can  be  evaluated  by  utilizing  (55)  to  yield  the  closed -form  solu¬ 
tion 


where 


P  =  exp(-a^)  ll-Q(bX,b)]  +^Q(a^c-r,  a^c+r) 
+  ^  (l  -  ^)  exp(-a^c)  , 


a  = 


r  2 +/?(!+ X*^) 

,  b  =  -===  ,  c  - - 


1-X 


/  2  2 


(98) 


(99) 
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For  r  =  0,  (99)  reduces  to  (H^/r)/2,  which  checks  the  last  equation  in 
Ref.  9.  For  M>2,  integrals  of  powers  Q  ar'  .equired;  in  this  case, 
generalizations  of  (59)  are  necessary. 


19/20 
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Preceding  page  blank 
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(8) :  Differentiate  (5)  with  respect  to  a,  and  use  (7)  and  (A -4). 

(9) :  Substitute  x  =  pt  in  (1)  and  reidentify  a  as  a/p,  and  b  as  bp. 

(10) :  Use  the  relation 


b  a>  (» 

dx  f(x)  +  I  dx  f(x)  =  J  dx  f(x) 


u  ao 

J  dx  f(x)  +  j 
0  b 


(A-11) 


and  employ  (9)  and  (2). 

(11):  Notice  first  from  (1)  that  Q(a,b)  is  an  analytic  function  of  a  and 
b  for  all  finite  a,b.  Then  replace  p  by  ip  in  (10),  and  utilize 


Q(iia,+ib)  -  Q(ia,ib),  a,b  real ,  (A-12) 

which  follows  directly  from  (1)  or  (3).  Also  (A-12)  is  real,  as  may  be  seen 
from  (3). 

(12) :  Substitute  a  =  bp^  in  (11)  and  use  (6). 

(13) :  Integrate  (9)  by  parts,  with 

2  2 

u  =  iQ(ax),  dv  =  dx  X  exp(-p  x  /2) ,  (A-13) 


and  then  employ  (5). 

(14) :  Use  (A-11),  (13),  and  (2). 

(15) :  Replace  p  by  ip  in  (14)  and  utilize  (A-12). 

(16) :  Put  a  =  bp^  in  (15)  and  use  (6). 

(17) :  Take  the  derivative  of  (9)  with  respect  to  a,  and  employ  (8). 

(18) :  Use  (A-11),  (17),  and  (2).  An  extension  of  (18)  is  also  supplied  by 
replacing  p  by  ip  (i.e.,  p^  by  -p2);  the  special  case  of  a  =  bp^  is  very 
simple,  upon  use  of  (6). 
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(19) ;  Take  the  derivative  of  (9)  with  reepeet  to  p2,  and  employ  (7)  aad 

(20) :  Use  (A-11),  (19),  (2). 

(21) :  From  Ref.  13  .  6.633  2,  we  have 


/  2  2\  <u 

b  exp^  -  J  ^  exp(-t' 


/2)  Jo(at)  J^(bt) 


the  last  step  by  Ref.  14,  9.1.30.  Integrating  (A -14) 
using  (10),  yields  ^  ^ 


_  d  r  2 

ab  j  dt  exp(-t  /2)  j^(at)b  J^(bt) ,  (A -14) 
with  respect  to  b,  and 


K  r  2  P 

b  J  dt  exp(.t  /2)  j^(at)  J  (bt)  =  f  dx  x  e: 

0  J 


2  2\ 

expl.^]l^(ax) 


^  ^  ~Q(a,b). 


(A-15) 


Let  t  =  px  in  (A-15)  and  reidentify  a  as  a/p.  and  b  as  b/p. 

(22) :  Replace  a  by  ia.  and  b  by  ib  in  (21).  and  use  (Ref.  14.  9.6.3] 

Jj(ix)=il^(x). 

(23) .  Put  b  =  a  in  (22)  and  omploy  (6). 


(M) 

j  dt  exp(-t)  l^fbt)  =  I  d,  exp,-t,  i  J  oos 


lO 


■I 
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_  exp  (-a)  r  exp(ab  cos  0) 

~  V  J  1  -  b  cose 
0 

r*  00  n  1 

_  gxR^-a),  I  exp(ab  cos0)  ^  jl  +  2  (“^)  cos(nfl)| 

0  L  1  J 

[00  n  n 

v->-E(v)vH 

n=l  -* 

=  “  [2Q(v/5:(r5),  ^(I^)  -  exp(-a)  yab)J  ,  (A-17) 

where  (4)  has  been  utilized  to  sum  the  series.  When  we  let  t  =  px  in  (A-17) 
and  reidentify  a  as  ap,  and  b  as  b/p,  (24)  follows. 

(25):  Use  (A-11),  (24),  and  (2). 

(26):  This  is  a  limiting  case  of  (25)  as  b — ►p-.  However,  (25)  ap¬ 
proaches  0/0  as  b— ►p-  (see  (6)).  Therefore,  let  b  =  Vp^-s2  in  (25)  and 
apply  L'Hospital’s  rule  at  s  =  0.  Upon  use  of  (7)  and  (8),  (26)  follows. 

(27):  Integrate  by  parts,  with 

u  =  exp(-px),  dv  =  dx  Ij^(bx) ,  (A-18) 


and  then  employ  (24). 

(28) :  Use(A-ll),  (27),  and  (2). 

(29) :  This  is  a  limiting  case  of  (28)  as  b — ►p-.  Apply  L' Hospital's 
rule,  as  in  the  derivation  for  (26),  above,  and  use  (7)  and  (8). 

(30) :  Take  the  derivative  of  (24)  with  respect  to  p,  and  utilize  (7)  and 

(8). 

(31) :  Use  (A- 11),  (30),  and  (2). 
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(32):  Define 


a 

Jdxx 


exp(-px)  I^(px).  n=0.1 


In  (A-19),  for  n  -  0,  integrate  by  parts,  with 


u  -  exp(-px),  dv  =  dx  X  I^(px) . 


Then,  using  {A-9), 


+A^. 

AdditionaHy,  in  (A-19),  for  n=0,  let  Instead 


U  =  X,  dv  =  dx  exp(-px)  I^(px) . 


Then,  using  (2G), 


(A-19) 


(A-20) 


(A-21) 


(A -22) 


A„  =  a‘  exp(-pa)  1  I^,(pa)  ^  i^ipa))  -  ,a^,.a^)  .  ,a-23) 

Equations  (A-21)  and  (A-23)  can  now  be  solved  for  both  A„  and  Aj. 

(33) :  Take  the  dorlvativo  of  (27)  with  respect  to  p,  and  utilize  (7)  and 

(34) -  Use  (A-11),  (33),  and  (2). 

(35) :  See  the  derivation,  above,  for  (32). 

(3G):  Integrate  by  parts,  with 

u  =  Q(ax,b),  dv  =  dx  X  cxp(-p“x^/2) ,  ^^-24) 

and  then  employ  (8)  and  (13). 

(37):  Use  (A-11),  (36),  and  (2). 


(38) :  Replace  p2  by  -p2  in  (37). 

(39) :  Integrate  by  parts,  with 

2  2 

u  =  Q(px,b),  dv  =  dx  X  expft)  x  /2) , 
and  then  employ  (8),  (2),  (A -4),  and  (5). 

(40) :;  Integrate  by  parts,  with 

u  =  Q(b, ax),  dv  =  dx  X  exp(-p\^/^) , 
and  then  employ  (7)  and  (9). 

(41) :  Use  (A-11),  (40),  and  (2). 

(42) :  Replace  p2  by  -p2  in  (40).  Integral  (42).  converges  for 
cause  as  b  —  m  ,  from  (1), 

<»  /  2  2\  ir~  / 

b  b 

=  J^4>(a-b)  ~  (2Tab)"^/^  exp^- 

wliere 

t 

4>(t)  =  J  dx(27r)  exp(-x^/2) . 

(43) :  Replace  p2  by  -p2  in  (41). 

(44) :  Integrate  by  parts,  with 

2  2 

u  =  Q(b,px),  dv  =  dx  X  exp(p  x  /2) , 
and  then  employ  (7),  (2),  and  (A-9). 
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(A -2  5) 


(A -2  6) 


p  <  a  be- 


(A-27) 


(A -2  8) 


(A -2  9) 
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(45):  Let  f  denote  the  left  side  of  (45).  Then,  using  (7)  and  Ref.  13, 
6.633  2,  there  follows 


/  2^2  2\ 
l-p  b  +c  \, 
■exp  — 5— ~  I 


(A-30) 


Since  f  -  0  at  b=  from  (45)  and  (2),  we  have 


T  /  2  2  2\ 

p  +  a  \2(p  +a  )/ 


and  (45)  follows  upon  use  of  (9). 


/  acx  \ 

V^a^j 


(A-31) 


(46) :  Employ  (5)  and  then  (9);  (2);  Ref.  13,  6.633  2;  and  (45)  to  evaluate 
the  resultant  integrals.  Next,  use  (5)  again. 

(47) :  Express  Q  in  integral  form  via  (1),  interchange  integrals,  and 
utilize  (A -9)  to  obtain 


no  /  2  2\ 

c 

“  c  ^  [I*  •  “p(-  ’(/v)] 


(A -32) 


where  we  have  employed  Ref.  IS,  6.633  2.  Equation  (47)  follows  directly  from 
(A-32). 

(48):  Integrate  by  parts,  with 

2  2  1  /b^  \ 

u  =  Q(c,ax),  dv  =  dx  exp(-p  x  /2)  Ij(bx),  v  =  -  expf— j  Q(px.  b/p) , 

'2n  / 


(A -33) 


and  employ  (7),  (8),  (2),  and  (45). 
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(49) :  Use  (5)  and  then  (13),  (2),  (22),  and  (48)  to  evaluate  the  resultant 
integrals.  Then  utilize  (5)  again. 

(50) :  Let  b  =  ac  in  (49)  and  use  (6). 

(51) :  Express  Q  in  integral  form  via  (1),  interchange  integrals ,  and 
utilize  (A -4)  and  Ref.  13,  6.633  2. 

(52) :  Take  the  derivative  of  (45)  with  respect  to  c,  and  use  (A -4)  and 

(8). 

(53) :  Take  the  derivative  of  (46)  with  respect  to  c,  and  employ  (7)  and 

(A-4). 

(54) :  Take  the  derivative  of  (47)  with  respect  to  c,  and  use  (A-4)  and 

(A-9). 

(55) :  Integrate  by  parts,  with 

2  2 

u  =  Q(ax,bx),  dv  =  dx  X  exp(-p  x  /2) ,  (A -34) 

and  employ  (2),  (7),  (8),  (24),  and  (27). 

(56) :  Use  (A-11),  (55),  and  (2). 

(57) :  Integrate  by  parts,  with 


u  =  Q(cx, bx),  dv  =  dx  X  , 


(A -35) 


and  employ  (7),  (8),  and  (2).  The  integral  Jv  du  can  be  reduced  to  the 
forms  of  (31)  and  (34)  by  the  substitution  y  =  x2/2 . 

(58) :  Use  (G),  let  y-x^,  and  utilize  (26). 

(59) :  Express  one  of  the  Q-functions  in  integral  form  via  (1),  interchange 
integrals,  and  use  (45)  to  evaluate  the  inner  integral.  Next,  use  (36)  to  evalu¬ 
ate  the  remaining  integral. 

(60) :  Express  Q  in  integral  form  via  (1),  and  interch:uige  integrals. 
After  taking  a  derivative  with  respect  to  a  ,  utilize  Ref.  13,  6.618  4. 
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6.618  ‘  ‘  ^  ““  ««f-  W. 

.-h.  E”P5oy  (13)  with  p=l.  to  express  Q  la  integral  lortn- inter¬ 
change  mtegrals;  use  Rel.  13,  6. 631  7;  and  then  (24)  and  (2^ 

(63) :  Use  (5);  Ref.  13,  6.618  4;  and  (62). 

(64) :  This  is  a  definition  of  4>. 

(65) ;  Integrate  by  parts ,  with 

u=4>(ax),  dv  =  dx  X  exp(-pV/2)  I^(bx) ,  (A -36) 

and  employ  (64),  (9),  (2),  and  (62). 

(66) :  Integrate  by  parts,  with 


u=«t.(-ax),  dv  =  dx  X  I^lbx) , 


(A -37) 


and^employ  (A-9)  and  ReJ.  13,  6. 618  4,  after  taking  a  derivative  with  respect 


(67):  Split  the  inte'yral  into  a  sum  of  integrals  over 
Let  y  -  -X  in  the  integral  over  (-«»,  0)  and  use 


(-'«,  0)  and 


(0,«>). 


4>('ax)  =  1  -  (f.(ax) , 

which  is  deducihle  from  (64).  Next,  employ  (9),  (2).  and  (65). 

(68) :  Integrate  by  parts,  with 

u  =  <|>(ax),  dv  =  dx  exp(-p^xV2)  Ij(bx) , 
and  employ  (64),  (13),  (2),  and  (63). 

(69) :;  Integrate  by  parts,  with 


(A -38) 


(A -39) 


u=<l>(-ax),  dv=dxlj(bx), 
and  employ  (A-4)  and  .Hcf.  13,  6.618  4. 


(A -40) 
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(70) :  Use  a  procedure  similar  to  that  above  for  (66),  and  then  employ 
(13),  (2),  and  (69). 

(71) :  Let  x  =  ^Jy  and  employ  Ref.  13,  6.644. 

(72) :  Use  Ref.  13,  6.6317,  with  v  =  1,  to  eliminate  Ii(cx2).  The  re¬ 
sultant  integral  involving  lo(cx^)  follows  from  (71).  Interchange  integrals  in 
the  remaining  double  integral,  and  use  (21)  and  then  (62). 


(73):  This  integral  is  a  limit  of  (72)  as  q— *-c+.  As  q-^c+,  the  right 
side  of  (72)  approaches  ■  ' 

where  t  =  .^q-c.  However,  by  an  approach  similar  to  that  given  in  (A -27),  we 
find  that 


Q(a,b)~  /- 4>(a-b)  as  a,b— 
*  3. 


(A -42) 


but  b  -  a  »  1  is  not  required.  Equation  (73)  follows  upon  use  of  Ref.  14, 

9.7.1. 

(74):  Define 


Then 


jr 

g  =  *  1  cbc  cos(nx)  expf-  r — r-^ - ). 

'’n’Tj  *^\l-b  cos  X/ 


(A-43) 


(A -44) 


using  Ref.  16,  vol.  2,  p.  81,  Eq.  10.  From  (A-43),  it  is  seen  that  gQ  =  1  at 
a=  0.  When  this  fact  and  (25)  are  utilized,  (74)  follows.  i 

(75):  Take  the  limU  of  (74)  as  b — ►!-  and  use  (A-42). 
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(76):  Define 


00 

C  2  2 

f^'=  I  dx  exp(-qx  )  Jj^(bx)  I^(cx  ) . 


(A-45) 


Express  Ijj  in  integral  form  by  use  of  Ref.  14,  9.6.19;  interchange  integrals; 
and  use  Ref.  15,  1.701^  to  obtain 


f^  =  ^  J  do  cos(n0)  1^1  -  exp^- 


q  -  c  cos  e/ 


(A -46) 


Next,  utilize  (74). 

(77) :  Take  the  limit  of  (76)  as  q — ►  C;+  and  use  (A -42)  and  (A -38). 

(78) :  From  (A-43)  and  Ref.  16,  vol.  2,  p.  81,  Eq.  10, 


f, 

da  bx  I  * 


r  .  1  -  b  cos  x-1  /  a 

j  1  -  b  cos  X  1  -  b  cos  X 


(A-47) 


From  (A-43),  =g0  at  a=0.  Then,  to  find  from  (A-47),  we  must  be 

able  to  evaluate  fdx  gQ,  which,  from  (A-43)  and  (74),  becomes 

'  'o' 

a  a  r 

I  dx  g^  =  Jdx  j2Q(av^,  -  exp^-  — 

i  0  0  * 

where  a  =  yjl-u /u,  P  =  y/l-Hi/n,  and  u  =  vi:?.  When  we  let  =  t  in 
the  first  integral  of  (A -48)  and  employ  (57)  and  (25),  (78)  follows. 


(79):  Take  the  limit  of  (78)  as  b—  1-,  and  use  (A-42)  and  Ref.  14. 


9.7.1. 


See  the  derivation  above  for  (76),  set  n  -  1,  and  employ  (78). 
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(81) :  Take  the  limit  of  (80)  as  q — ►  c+,  and  use  (A -42)  and  Ref.  14, 
9.7.1. 

(82) :.  Express  Q  in  integral  form  via  (1),  interchange  integrals,  and 
use  (19)  and  (1). 

(83) :  Use  (A-11),  (82),  and  (2). 

(84) :,  Use  (5),  (83),  and  (1). 

(85) :  Use  (5),  (82),  and  (1). 

(87) :  Litegrate  by  parts  on  (86),  with 

/  2  2\ 

m+n-1  /  X  +a  \  ,  ,  -n+1 ,  .  ,  ,. 

u  =  X  expl - - — •),  dv  =  dx  X  *  (A -49) 

and  employ  (A -4)  and  (86). 

(88) :  Integrate  by  parts  on  (86),  with 

/  2  2v 

u  =  x"^  I^(ax),  dv  =  dx  X  expf-  ^  ^  ^  j,  (A -59) 

and  employ  (A-y)  and  (86). 

(91) :;  Express  Q  in  integral  form  via  (1),  interchange  integrals,  and 
use  (86)  and  (9). 

(92) ;  Express  Q  in  inlegral  form  via  (13),  with  p=l,  inteichange 
integrals,  and  use  (86)  and  (13). 

(93) :  Integrate  by  parts,  with 

u  =  <p(-ax),  dv  =  dx  Ij^(bx; ,  (A-51) 

and  employ  (A-4)  and  (86). 
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